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LOCAL SMOOTHING ESTIMATES FOR THE MASSLESS 
DIRAC-COULOMB EQUATION IN 2 AND 3 DIMENSIONS 

FEDERICO CACCIAFESTA AND ERIC SERE 


Abstract. We prove local smoothing estimates for the massless Dirac equation 
with a Coulomb potential in 2 and 3 dimensions. Our strategy is inspired by [9] 
and relies on partial wave subspaces decomposition and spectral analysis of the 
Dirac-Coulomb operator. 


1. Introduction and generalities 


The massless Dirac equation with an electric Coulomb potential reads 


( 1 . 1 ) 


iu t + V n u — —rU = 0, u(t, x) : M f x 
\x\ 

u{0,x) = f(x) 


-> C 


N 


where the massless Dirac operator T> n is defined in terms of the Pauli matrices 

d.2) *■ = (! i). ^=(°i o')’ ff3= (o -l) 

as 

T>2 = —i(a!d x + a 2 d y ) = ^ o ^ ) 

(we denote d z = d x — id y and d~ = d x + id y ) in dimension n — 2 with N = 2, and 

3 

P.3 = VI ai k d k = —i(a ■ V) 

k =1 

where the 4x4 Dirac matrices are given by 

(1.3) ^ k — 1,2,3 

in dimension n = 3, with N = 4. 

The Dirac equation is widely used in physics to describe relativistic particles 
of spin 1/2 (see e.g. [29], [40]). In particular, the 2D massless equation is used 
as a model for charge carriers in graphene, a layer of carbon atoms arranged in a 
honeycomb lattice (see e.g. the physics survey [33] and the mathematical papers 
[22], [25], [30]). Note that there have been many rigorous results in the recent 
years on the stationary (massive) Dirac equation, essentially in dimension 3 (see 
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e.g. [18] for references). Comparatively, its dynamics has been less investigated 
up to now. 

The aj and oij matrices were introduced in view of making the Dirac operator 
a square root of the Laplace operator: therefore, they satisfy by construction the 
following anticommutating relations 


T , j, k 1, 2, 

otjC%k T a k a.j 25^14, j, A 1,2,3, 

(we recall that the spectrum of T> n is the whole line M). These conditions ensure 
that 

(id t - V n )(id t + V n ) = (A - d^ t )I N , 

which strictly links the dynamics of the free massless Dirac equation to a sys¬ 
tem of N decoupled wave equations. Therefore, dynamical properties for the free 
Dirac equation can directly be derived from their wave counterparts: dispersive 
properties of the flow can be mainly encoded in the celebrated family of Strichartz 
estimates , which are given by 


(1.4) 



< II/IIl* 


where the exponents (p, q ) are wave admissible, i.e. satisfy 


2 d — 1 

- H- 

P Q 


d- 1 
2 


2 < p < oo, 


2{d — 1) 
d -3 


> q > 2 


and e lWn represents the propagator for the solutions to (1.1) with V = 0 for 
n = 2, 3. We stress the fact that the estimate corresponding to the couple (p, q) = 
(2, oo) in dimension 3, the so called endpoint, that is 


(1.5) 


a itV 3 


/II 


i4l% 


< 




is known to fail as the corresponding one for the 3D wave flow. In [31] a refined 
version of this estimate involving angular spaces was given. As a consequence, the 
well-posedness of the celebrated cubic nonlinear Dirac equation was established, 
assuming initial data in H 1 with slight additional regularity on the angular vari¬ 
able. The problem was afterwards completely solved in [4], See also [5] for the 2D 
case. 

Another family of a priori estimates encoding informations about dispersion is 
given by the so called local-smoothing type estimates, which turn to be particularly 
useful to handle potential-type perturbations. Different versions of these estimates 
are available: to the best of our knowledge the most general one in this setting is 

(1.6) IK 1/2 e i ®>/IU; i | < C||/||i* 


where w a (x) = |x|(l + |log|x||) <r , a > 1. Both (1.4) and (1.6) can be found in 
[14] (see also [17], [7]). We are not aware of any corresponding result for the 2D 
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case; nevertheless, the strategy developed in [14] can be suitably adapted to this 
context to obtain similar results. 

In recent years a lot of effort has been spent in order to investigate what happens 
to dispersive flows when potential perturbations come into play. In [14] both (1.4) 
and (1.6) were proved for the flow e lt ^ Di+x 1 under the assumption 

\V(x)\ <— cr>l, n = 3, 
w a {x) 

for some 6 sufficiently small, which by the way is a sufficient condition to guarantee 
that the operator T>^ + V is selfadjoint (see [40]). Similar results are available for 
magnetic potentials as well (see e.g. [8], [13]). The angular endpoint estimates 
proved in [31] have been generalized in [11], [12] to small potential perturbations, 
by introducing some new mixed Strichartz-smoothing estimates. At the same time, 
some effort has been spent in order to find examples of potentials such that the 
correspoding flows do not disperse (in the sense discussed above). This problem 
has been tackled for the magnetic Dirac equation in [2] in 3D, essentially showing 
that for vector potentials of the form A{x) = \x\~ s Mx with 6 G (1,2) most of 
the mass of the solution is localized around a non-dispersive function. These 
arguments suggest that, heuristically, the degree of homogeneity of the operator 
works as a threshold for the validity of dispersive estimates: therefore, for the Dirac 
equation, the Coulomb potential is to be thought of as a critical case (we stress the 
fact that the restriction to the massless case seems unavoidable now). In the very 
last years the problem of understanding dispersive estimates for scaling invariant 
potentials has been approached for other dynamics: in [9] and [10] the authors 
have proved indeed that Morawetz and Strichartz estimates can be recovered for 
both the Schrodinger and wave equations with inverse square potentials, and in 
[21] the L 1 — y L°° time decay estimate is proved in the electromagnetic case. 

The aim of this paper is to adapt the techniques of [9] to the massless Dirac 
Coulomb equation, and thus to prove local smoothing estimates for solutions to 
(1.1). The strategy of proof of [9] is fairly straightforward: the basic idea is to use 
a spherical harmonics decomposition to reduce the equation to a radial one and 
then rely on the Hankel transform to diagonalize it. Several problems arise when 
trying to adapt these ideas to the Dirac setting. The first difficulty is given by the 
fact that the Dirac operator does not preserve radiality, and therefore one needs 
a more complex, but still classical, spherical harmonics decomposition related to 
the SU 2 group action (see [40], section 4.6). Another issue is that we cannot use 
the Hankel transform, and we have to build up a relativistic analogue. To do this, 
we need to study the continuous spectrum of the Dirac-Coulomb model, both in 
2 and 3 dimensions. The corresponding generalized eigenfunctions are well known 
(see for instance [29] and [34]), but their rich structure will raise some additional 
technical difficulties. We stress the fact that, to the best of our knowledge, this 
is the first result concerning dispersive properties of Dirac-Coulomb, an operator 
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which plays a fundamental role in relativistic quantum chemistry. We hope that 
our result will be useful in the study of the dynamics of nonlinear models involving, 
for instance, several relativistic electrons in a molecule. 

Before stating our main result, we introduce some notations that will be used 
throughout the paper. 

Notations. 

We shall use the standard notation H s for the homogeneous Sobolev space with 
the norm \\f\\fj B = || |-D | S /|| L 2 where \D\ = (—A) 1 / 2 , and L p t L q x = L p (R t ; L q (W£)) 
for the mixed space-time Strichartz spaces. We denote with 12 s the operator 

(Ll s 4>)(x) = \x\ s (f){x) 

and with a little abuse of notation we will use the same symbol to indicate the 
operators which are pointwise equal for all times, 

(Ll s fj)(t, x) = \x\ s i/>{t, x). 

We introduce the following sets of indices 

(1.7) B 2 (k) := {k : k = 1/2 + m, m E Z} 

(1.8) B 3 (k) :={k:kE Z\{0}}. 

As most of the forthcoming objects will be significantly different in 2D or 3D, 
we will often use the apex n to distinguish the dimensions; this will not have to 
be confused with standard powers, but interpretation will be clear from time to 
time. For the definition of the angular spaces Ti/, n = 2, 3 we refer to forthcoming 
Proposition 2.1 and subsequent Remark 2.1. 


We are now ready to state our main Theorem. 


Theorem 1.1. Let n = 2,3, v n E (—and let k n E B n (k). Let u be a 
solution of (1.1). Then, for any 


1/2 < a < ^^2 + 1/2 

and any f E L 2 ((0, oo)) <E> TLf kn there exists a constant C = C{v,a,k) such that 
the following estimate holds 


(1.9) 



^ Di 

T> _ — 

^ n | | 

1/2-a 

cr a 

U 


\ X \ 



D t Ll 


<C\\f\\ L 2. 


Remark 1.1. The dependence on the dimension of u, which will be neglected in 
the rest of the paper, has been added in the statement of the Theorem to allow us 
to remark the different ranges of admissible charges in dimensions 2 and 3. These 
ranges are in fact imposed by the condition on the Dirac Coulomb operator to be 
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self-adjoint, which is indeed guaranteed in those intervals (see Remark 1.4). We 
stress also the fact that for v approaching the endpoints, respectively ±1 in 3D 
and ±1/2 in 2 D, the range of admissible exponents for the lowest value of k n in 
estimate (1.9) shrinks to zero. 

Remark 1.2. It seems to be possible to generalize our proof of Theorem 1.1 to any 
space dimension n > 2. As it will be clear, the key ingredients are given indeed 
by the construction of a suitable integral operator which transforms the prob¬ 
lem into an ODE, and then by the careful analysis of some interaction integrals 
between generalized eigenstates of the Dirac-Coulomb equation that naturally ap¬ 
pear. Therefore, the crucial step consists in finding the generalized eigenstates of 
the Dirac-Coulomb equation: this has been performed in [24] in higher space di¬ 
mension n > 4 by generalizing the spherical wave decomposition (see Proposition 
2.1), i.e. by carefully analyzing the symmetries of the SO(n ) group. Anyway, we 
prefer not to deal with this problem here both because technicalities would lead 
us too far and because of the scarce physical interest of it. 

Remark 1.3. The endpoint case a = 1/2 in (1.9) would permit to recover with 
a standard argument the full set of Strichartz estimates for the perturbed flow 
(see e.g. [9]); even though we don’t have a concrete counterexample, this estimate 
seems to fail. Forthcoming remark 3.1 gives more details on the issue. 

Remark 1.4. In order to have a unitary flow for equation (1.1) one needs the 
operator D n — py to be selfadjoint; many papers have been devoted to the study 
of this property (see e.g. [37], [42], [32], [27], [19], [20], [3]) in the 3D setting (note 
that the methods introduced in these papers are quite general and can be adapted 
to other dimensions). In 3D the essential self-adjointness is in fact guaranteed 
if \v\ < while in the interval ^ < \v\ < 1 it is still possible to build a 
distinguished selfadjoint extension. In the 2D case, as soon as v is nonzero, the 
Dirac-Coulomb operator V- 2 — defined on (^(M 2 ) is not essentially selfadjoint 
(this is an immediate consequence of Theorem 4.16 in [40]). But C. Warmt ([41], 
Satz 2.2.6) recently proved the existence of a distinguished self-adjoint extension 
if \v\ < 1/2. Therefore, in the range of u n for which we prove our Theorem 1.1, 
the operator V n — py is self-adjoint (in the sense just discussed). 

The plan of the paper is the following: in section 2 we build the setup, review¬ 
ing the theory of partial wave subspaces and the spectrum of the Dirac-Coulomb 
operator in 2D and 3D; then we define our analogue of the Hankel transform used 
in [9]. Section 3 is devoted to the proof of Theorem 1.1. A short appendix contain¬ 
ing some generalities on special functions is included at the end for the reader’s 
convenience. 

Acknowledgments. The first named author is supported by the FIRB 2012 
” Dispersive dynamics, Fourier analysis and variational methods” funded by MIUR 
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2. The setup. 

We devote this section to recall some classical facts about the radial Dirac oper¬ 
ator and his spectrum, including a short discussion on how to explicitly calculate 
the generalized eigenstates of the continuous spectrum. For the sake of brevity we 
will omit several details, referring to [29] and [40] for a complete discussion and 
for more informations on the topic. Also, we construct an integral transform that 
will be crucial in our proof. 

2.1. Partial wave decomposition and continuous spectrum. The main in¬ 
gredient we need to introduce is the so called partial wave decomposition. We 
collect in Proposition 2.1 the tools we will need in what follows that is, essentially, 
the fact that the Dirac-Coulomb operator can be seen as a radial operator with re¬ 
spect to some suitable decomposition, both in dimensions 2 and 3. We will denote 
for brevity with A 2 and .A 3 the following families of indices 

(2.1) A 2 1/c G Z + - j> ; 

(2.2) A 3 := ■ ■ ■, rrij = -j,-j + 1, ..., +j, k 3 = ±(j + 1/2)1 

Moreover, we define the 4-vectors 


where 

a ™, j_ / yi + yjjifi 2 ) 

J - 1/2 vv ( ) 

0 m j 1 ( \/j + 1 - W/F ) 

3+1/2 V2J+2 ( + 1 - m, J 

and Y, m are the usual spherical harmonics. 

Proposition 2.1. Let n = 2,3. Then, we can define unitary isomorphisms be¬ 
tween the Hilbert spaces 

L 2 (M 2 ) 2 = T 2 (( 0 , 00), dr) (^) h\, 

k£A2 

L 2 (R 3 ) 4 = L 2 ((0,oo),dr) (g 

j,mj,kj&A3 


""U.TO'T 1 / 2 ) 


j±l/2’ 


(2.3) 


7 + 

~‘m j ,^(j+ !/ 2 ) 


itt 


I'T 1/2 

0 
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by the following decompositions 


(2.4) 


*»(*> = £ ^ 
keA 2 v 


F k [r)e i{k ~ l/2)e 

Gk(r)e l ^ k+1 ^ e 


) 


(2.5) 


*3(1) = X] 

j,m.j,kj£A3 


2 \/ 7T V ( r )^j,mj,kj (^> ^0 


V )^j,mj,k 


which hold for any <3> n G L 2 (M ra ,C iV ). Moreover, the Dirac-Coulomb operator 
leaves invariant the partial wave subspaces C^°((0, 00 )) ® hf and, with respect to 


the basis |e^ fc V 2 ) 0 , e*( fe+1 / 2 ) 0 } an d 
by the radial matrices 


-‘j,rrij,kj'> ' 'j,rrij,kj 


is respectively represented 


( 2 . 6 ) 4 = 


A. _i_ h 

dr r 


d 1 /c 
dr r 


d 3 

, m j ,kj 


_d_ 

dr 


1—ki 


d I 1+fcj 
dr r 


TTie Dirac-Coulomb opeator V n — - on Cfi°(M. n ) N is unitary equivalent to the 
direct sum of ’’partial wave Dirac operators” d, hm j.., 

7* ©4. *-7 


®2-7-©4. ®3 - - = © <f 

For /rred k 2 G B(k 2 ) and k 3 G B{kf) we denote with 

(2T «k= © 

k: | fc| >| fc 2 1 i 0 < I &31 -1/2,mj ,kj: 

Proof. See [40]. 


3 

j,m k ,k k 


hj, mj ,kj- 


□ 


Remark 2.1. Notations in Proposition 2.1 are very heavy; this is partly due to 
the richness of the Dirac operator itself and partly to the fact that the spherical 
coordinates decomposition changes in dimensions 2 and 3. We should stress the 
fact that formulas (2.6) only depend on k: therefore, in what follows we will abuse 
notations neglecting the dependence on the other parameters in the case n = 3 in 
order to have a unified presentation. 


We now show how the functions of the continuous spectrum for the Dirac- 
Coulomb operator are obtained. We want to deal with the eigenvalue equation 

^ n = e^ n , e > 0, n = 2,3. 

The application of Proposition 2.1 suggests to seek for solutions in the form, for a 
fixed value of e, 
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( 2 . 10 ) 


t,3 / ( Fl (r)Ef(x) 

4 7tJr,x) = Jfe e \ 

’ V G lA r )-k( x ) 


n = 3 


so that equation (2.8) will in fact reduce to a sytem of ODE for the radial compo¬ 
nents Fjf(r), G%(r) which is the same in both dimensions, apart for some constants, 
and can be expliticly solved (for the details on these calculations we refer to [29] 
section IV for the 3 D and to [34] for the 2D case). We thus have the following 
formulas, for fixed values of k and e > 0, and n = 2,3, 

( 2 . 11 ) 

+ 1 + “d{e*< Er+ « 1 F 1 (7 - ii/,27 + 1, -2 ier)} 


T?n _ 

** k.e 


r(2 7 +1) 

l ^^T{2++l) ^ eW2 ( 2gr ) 7 ~ L ^ Im {e* (£r_K) iFi( 7 - w, 2 7 + 1, -2 ier)} 

where i Fi(a,b,z) are confluent hypergeometric functions (see Appendix for the 
definition and some properties), 7 = \Jk 2 — u 2 and e~ 2t ^ is a phase shift. 

We will denote with 


( 2 . 12 ) 


G n 

k,e ~ 


(2.13) 


.„ . , , FT e (r) 

^W=( Gi < r ) )• 


n = 2, 3 


the vector of radial coordinates of the generalized eigenfunctions. 

Notice that with this choice the following normalization condition holds for 
(2.9)-(2.10) 


(2.14) 


'•A 00 _ — 

'if n flr)f; 1 f,(r)r n - 1 dr = 


c n— 1 


n = 2,3 


meaning with that that the function 

K(r,r') = r”^ 1 j e^ le fj n e (r)'ipf,(r)£ n ~ 1 d£ 
is a ^-approximation for l —y 0. 


Remark 2.2. We point out that, due to our choice of normalization (2.14), the 
homogeneity of the functions -0 n (r) with respect to £ and r is the same; this fact, 
that is not true anymore in the massive case, will be crucial in what follows, and 
therefore prevents the application of our strategy to the massive case. 


Remark 2.3. The sign of the Coulomb potential plays in fact a virtually unessential 
role in what is above, and the wave functions for repulsive fields can be obtained 
simply by changing the sign of v in (2.11)-(2.12). 

Remark 2.4. Formulas (2.11)-(2.12) give the ’’eigenstates” of the continuous spec¬ 
trum corresponding to positive energies £ > 0. Using a charge conjugation argu¬ 
ment (or directly working on the explicit radial equations), it is possible to obtain 
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the corresponding formulas for the radial components (2.9) for negative energies: 
it turns indeed out that the function rj)_ k ^ e (r) solves 

(2.15) d k if_ Ke {r) = n = 2,3, 

where the operator d k is the one defined in (2.6) with opposite sign of the charge v. 
This shows that the ’’eigenstates” corresponding to negative energies for equation 
(2.8) can be understood by mean of (2.15) to be ’’eigenstates” corresponding to 
positive energies, with opposite value of j and opposite charge. Therefore, denoting 
with F, G the functions obtained by (2.11)-(2.12) by changing the sign of z/, we 
obtain for the radial coordinates of the ’’eigenstates” corresponding to negative 
energies as in (2.13) the following relation 




jfe, —£ 


r = 


F k,-e( r ) 


F -k,e( r ) \ 

G-k,e( r ) J ' 


2.2. The integral transform. We now introduce the crucial integral transform 
that will be used in the proof of the main result, consisting in a projection on the 
continuous spectrum of the Dirac-Coulomb operator. Throughout this subsection 
we will omit the dependance on the dimensions of the functions, as the calculations 
carried out will be the same. 


Definition 2.2. Let $ G L 2 ((0, oo), dr) ® h k for some fixed k and let ip(r) = 
(</?i(r), </?2 (t)) be the vector of its radial coordinates in decomposition (2.4)-(2.5). 
We define the following integral transform 


(2.16) 




F k M e ) 


/ 0 +oo ^fc, e (r)v?(r)r n 1 dr 

C (/o + 0 °^fc-e( r )^( r ) rri_lrfr ) 


r»+oo 


= / H k {sr) ■ Lp(r)r n v dr 

Jo 

where we have introduced the matrix (keep in mind Remark 2.4) 


(2.17) 


H k 


F k,e ( r ) G k ,s( r ) \ 
F k ,-e(r) G k - e (r) J 


We collect in the following proposition some important properties of the operator 

V k . 


Proposition 2.3. For n = 2,3 and any $ G L 2 ((0, oo), dr) <8) h k the following 
properties hold: 

(1) V k is an L 2 -isometry. 

(2) V k d k = a^FtV k - 

(3) The inverse transform ofV k is given by 



(2.18) 


Ht(er) ■ ¥>( e)e"~'de 
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where H* — ( FkA r ) Gk,-e{r) \ 
where H k -y Gfee(r) Fk _^ r) )■ 

(4) For every a E M. we can define the fractional operators 


(2.19) 


r*+oo 


A kfk(r) = V k a 3 Q a V k \ p k (r ) = / S k (r,s) ■ ip k (s)s n 1 ds 


where the integral kernel S k (r, s ) is the 2x2 matrix given by 


( 2 . 20 ) 


+oo 


S k (r, s ) = / H k (er) • H* k (es)e n ~ 1+a de 


Jo 

Remark 2.5. Property (1) allows, by standard arguments, to extend the definition 
of the operator V k to functions in L 2 . 


Remark 2.6. When summing on k, property (2.19) defines in a standard way 
fractional powers of the operator \T> n — j^|, which are used in the statement of 
Theorem (1.1). 

Proof. Property (1) and (2) come from the definition of V k , once noticed that 

V k (d k <p) = (lpk,±e, d k <p) = {dki>k,±eiF>) = ±e(^fc,± e,p) 

and from normalization relation (2.14). Property (3) is a consequence of (1) (HI 
is the adjoint of H k ). 

To prove property (4) we use the definition of A k to write 
A kPk(r) = V k FL a Vfi l y k (r) 

/ +oo / r+oo \ 

H k (er)e n ~ 1+a (j Hl(es)^ k (s)s n - l ds\ de. 
Exchanging the order of the integrals yields (4). □ 


3. Proof of Theorem 1.1 

We here prove our main result. Following [9], we rely on partial wave decom¬ 
position and on Proposition 2.3 to explicitly write the solutions on the single sub¬ 
spaces; then, a detailed and careful analysis of the kernel S k (r,s) of the operator 
A k defined in (2.19) will be performed to sum back. 

3.1. The strategy. We start with decomposition (2.4)-(2.5) and work on a single 
spherical space which will be identified by k G A 2 as in (2.1) and by the triple 
{j, m k , kj} G A 3 as in (2.2) in the case n = 3. With a slight abuse of notations, we 
will again neglect the dependence on everything but k in the 3 D case (as noticed 
in Remark 2.1, the radial Dirac-Coulomb operator in dimension 3 only depends 
on k = kj , see (2.6)), as well as the dependence on the dimension for the various 
objects, in order to perform much readable and unified calculations. We thus rely 
on Proposition 2.1 aiming to prove estimate (1.9) on a fixed partial wave subspace 
with a constant C bounded with respect to k, and eventually sum back. We thus 
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take an initial condition / with angular part in h k for some fixed k, and denote 
with L k f the solution to the initial value problem 


(3.1) 


iu t + d k u = 0, 
u(0, x) = f{x) 

with d k as in (2.6). We apply operator V k , which is an L 2 -isometry, to the LHS 
side of estimate (1.9) and use Proposition 2.3 to obtain (the application of the 
matrix 03 does not alter the L 2 norm) 

linw“|4l 1/2 -“n/|| iW = \\A^n'/ 2 -°v k L k f\\ LtLi 

where we have used (2.19). Now, VkL k f solves (see (2) of Proposition 2.1) 

id t V k L k f + a 3 nV k L k f = 0, 

V k L k f (0,£) = V k m, 

so that the solution to this problem is explicitly given by 

V k L k f {t,C) = e^V k m- 

Fourier transforming in time (since Tt-* T is an L 2 -isometry which commutes with 
A k ) gives (with a minor abuse, we will neglect the <73 matrix from now on as we 
will only need to pick the L 2 norm and thus the sign will be unessential) 

{J r tV k L k f)(T,Z) = (Vkf){p)8{T + Z). 

We can then write 

P+OO 


(3.2) 




s k “(f. S ) S ( T + s)V t f(s)s- 


~ds 


= S k a (Z,T)Vkf(T)T- 


I11 view of proving (1.9) we now need to calculate the Lr norm in time and space 
of the quantity above which then gives (notice that, as the angular part in decom¬ 
position (2.4)-(2.5) is L 2 -unitary, we only need to consider the radial integrals) 

P+OO /»+oo 


(3,3) 


((Ptf)’( T )St°(p,T) ) ■ (S k “(p,r)(P t f)(r))T 


2n—2a—l „2 


p dp. 


Since S%(p, r) 1 = S%(r, p), the integral in dp yields S k 2 q (t, t) and we are therefore 
left with 

p-1-00 


(3.4) 


wi'Mriviwiwr 


2n— 2a— 1 


dr. 


To conclude the proof of the estimate on the k-th space we need to bound 
(3.4) with the L 2 -norm of V k f (which we recall to be the L 2 norm of /). Using 
forthcoming Proposition 3.1 we thus have 

P+OO 


(3.5) 


(njY(T)S^(T,T)(V t f)(T)T 


2n—2a—l 


dr 




12 


FEDERICO CACCIAFESTA AND ERIC SERE 


r +oo 

< / Tr(^- 2 “(r,r))|(n-/)(r)| 2 r“- 1 dr 

J o 

/ +oo 

I(^/)(t)| 2 (t)t"- 1 dr 

= C k \\f\\ L 2. 

with C k a bounded constant of k. Expanding / in partial waves and using the tri¬ 
angle inequality together with the orthogonality of spherical harmonics concludes 
the proof. 


3.2. The main integral. The crucial step consists thus in evaluating 

r+oo 

F k , £ {ry + G k , e (rYe" 


Tr(5' -2a (r, r)) = Re 
(3.6) + Re 


\2 | ^n—l—2oL^^ 


IJ 0 


r»+oo 


F fc ,_ £ (r ) 2 + G fc ,_ £ (r) 2 e"- 1 - 2 “de 
U o 

= Re[/i(r) + .^(t)] . 


More than giving an explicit solution of such integrals, which seems to be possible 
but significantly complicated, we are more interested in proving the following 


Proposition 3.1. Let n = 2,3 k n € B(k ) as defined in (1.7)-(1.8) and u n as in 
Theorem 1.1. Then for every 

1/2 < a < - Vl + 1/2 

there exists a constant C kn such that 

Tr(Sj/ 2 “(r, r)) = C kl T-”+ 2 “ 

Moreover, sup \C kn \ < oo 


Proof. By recalling the structure of the generalized eigenstates (2.11)-(2.12), we 
have that integral (3.6) takes the form 

(3.7) h(r) = C k , u r 2 ^ n ~V 


x 


with 



e 2 7- 2 a e 2i6T iFi ( 7 _ 2 7 + 1 , -2fer) 1 F 1 (7 


iv, 2y + 1, —2 ier)de 


C k 




227+2 

7T 


|E(7+ 1 +H 2 2i£ 

r(2 7 + i ) 2 


where we recall that 7 = \Jk 2 — u 2 . As the only difference between I\ and I 2 is 
encoded in the sign of v (see Remark 2.4) which is virtually unessential, we will 
limit our discussion to I\, the other case being completely analogous. 

Integrals like (3.7) have already been object of study in literature (see e.g. [28], 
[23], [6]) since they are connected with some different physical problems, as for 
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instance matrix multipole elements for the Dirac Coulomb model, and they present 
some difficulties being not convergent at infinity as well as in the origin (see (4.3)). 
Singularity at the origin can be dealt with by a standard integration by parts 
argument (as the one used to define T function for negative real parts) and is 
discussed in [39] provided 2y — 2a + 1 is not zero or a negative integer, while 
singularity at infinity is slightly more difficult to be handled. To make things 
work, we need first of all to define Rfr, s) off-diagonal and then use a limiting 
procedure. Also, the introduction of a term of the form e~ Se with 5 E M + in the 
integral is needed in order to make it convergent; we will then denote with If the 
integral 

I s k (r,s) = C k ,„T 2 ^ n ~V 


x 


+oo 


^ 27 — 2a^e(i(r+s)—S) 


1 F 1 ('f - iv, 2y + 1, —2fer) 1 Fi(7 


iis, 2y + 1, —2 ies)de. 


Jo 

This modified integral can be explicitly computed, and its solution is written in 
terms of a double Appel series, which unfortunately turns out to be not convergent 
for the values we are interested in. To deal with this problem, it is necessary to 
introduce an analytic continuation of such a series: we don’t include the proofs of 
these results here as they can be found in corresponding references. 

We begin with the following 


Lemma 3.2. Let a, fi such that |a| + |/3| < \h\ and d not a negative integer nor 
zero. Then 

(3.8) f f d_1 e _,lt F(ai, £7, at)F(a 2 , b 2 , f3t)dt = h~ d T(d)F 2 (d; op, a 2 ; h, b 2l y, y) 
Jo h h 

where F 2 is defined in (4.5). 

Proof. See e.g. Lemma 1 in [36]. □ 

We should now point out that the ray of convergence of F(a; a 1; a 2 ; b\, b 2l x, y ) is 
|x| T \y\ < 1, which provides a first obstacle to the application of this result to our 
case. To overcome this difficulty, we rely on the following analytic continuation of 

f 2 . 


Lemma 3.3. If Re(b\ — cp) > 0, Refafi) > 0 
(3.9) F 2 (d] ai, a 2 , 61 , b 2 ] 1,1) 

^)WT + i) 


= e 


+e v 


r(a 1 )r(d — an +1) 
rfbfirfd-h + i) 


T(&i — afiTfd + Gp + 1 — b\) 
Proof. See [6]. 


3 F 2 fd, b 2 — a 2 , d — bi + 1] b 2l d — a\ + 1; 1) 
2 F 2 fd, a 2l d + 1 — b\] b 2l d + a + 1 — £7; 1) 


□ 
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For the sake of generality we have stated this Lemma with general values of the 
parameters; we list here for convenience the necessary substitutions to treat our 
case: 


(3.10) 


2 a; 


d = 2y + l 
ai = a 2 = 7 - iv\ 
bi = b 2 = 2y + 1. 

Since necessary conditions of Lemma 3.3 are fullhllcd with these values, this 
result allows us to pass to the limit l/(r, s ) for 5 —> 0, s —> r =: r, and thus leaves 
us with 

h { r ) = C fci! ,r 27 - (n " 1) (2ir)- (27+1 " 2Q) r(27 + 1 - 2a) 
xF 2 (2y + 1 - 2a; 7 - iv , 7 - iv , 27 + 1, 27 + 1; 1,1) 

with 


Ck,u = 


2 27+2 „|r ( 7 + i + H 




7r r(2 7 +1) 2 

Using now representation (3.9) for F 2 and plugging in values (3.10), we are led to 
(3.11) 
where 


h{r) = r 2a ~ n 

— e in{rt- 2 a-iv) 

E 


C k)V (2i)- { ~ w ~ 2a) Y(2 1 + 1-2 a){A + B ) 

r( 2 7 +l)r(l- 2 a) 


x 


r(7 — iv)T{j + 2 — 2a + iu) 
(27 + 1 — 2a) m (7 + 1 + iv)m{ 1 — 2a) m 


(27 + l) m (7 + 2 - 2a + iv)mm\ 


= e iTt(n-iv) 


r(2 7 + l)r(l-2a) 


x 




r(7 + 1 + iv)T{j + 1 — 2 a — iu) 
(27 + 1 - 2a) m (7 - iu) m ( 1 - 2a) m 

iv) m m\ 


m (27 + l)m( 7 + 1 — 2a 

Both the 3+2 functions above turn to be convergent: we can indeed explicitly 
check convergence condition (4.2) to give in both cases a > 0. Notice also that 
assumption a < 7 + 1/2 guarantees boundedness on the term T(27 + 1 — 2a). To 
conclude with, we now need to take the real part of (3.11) and show that it is a 
bounded function of 7. After some simplifications, we have 


Mh(r)) = 


2 l+2a e 2nu + i _ 2a)r(l - 2a) 


x Re 


7 T r(2 7 +1) 

^Z(2a + i) r(7 + 1 + iu )(7 - iv) (27 + 1 - 2a) m (7 + 1 + iv) m { 1 ~ 2a) r 
r ( 7 + 2-2 a-w) ^ (27 + l ) m (7 + 2 - 2 a + iv) m m\ 


+ 


e if(2a-l) 


r (7 + 1 — iv) 
r (7 + 1 — 2 a — iv) 


E 


(27 + 1 - 2 a ) m (7 - iv)m{ 1 ~ 2 a) r 
(27 + l ) m (7 + 1 - 2 a - iv) m m\ 
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where we have written j-( 2 "/+ l - 2a ) = e -*f (27+1-2+. Using now the basic property 
of the Pochhammer symbol (x) m = r(x + m)/T(x) many other terms simplify, 
leaving us with 


R e(h(r)) = 


21+20:^27 rv 


7T 


xRe 


e -if(2a+l) 


(7 - iv) Y 

77 

+e *f (2«— 


r(2y + 1 — 2 a + m)r (7 + 1 + m + iv)T(l — 2 a + m) 
r(2y + 1 + m)r (7 + 2 — 2 a + m — iv)m\ 

r(2 7 + 1 — 2 a + m)r (7 + m — iv)T{l — 2 a + m) 


r(2 7 + 1 + m)r (7 + 1 — 2 a + m + iv)m\ 
Stirling formula for Gamma function 

/27T f z\ z 


Y{z) = 



©■ph;)) 


shows that both the series above are asymptotic (in 7 ) to 7 1 , and therefore the 
proof is concluded. 

□ 


Remark 3.1. It should be noticed that when a —> 1 the term T(1 — 2a) produces a 
singularity of order 1 , as it is seen by applying well known property r(z)r(l — z) = 
n / sin( 7 rx). Anyway, it can be checked that this singularity vanishes when taking 
the real part. Calculations are very easy in the unperturbed case (i.e. when 
v = 0 ): in this case indeed everything is real except for the exponential terms, 

• 3 1 

which become respectively e _, 2 f and e l 2 7r which are purely imaginary. In the 
perturbed case things are a bit more complicated, but the problem has already 
been dealt with in literature by some limiting argument (see [ 6 ] and references 
therein). On the other hand, as a approaches 1/2 the singularity produced by 
the term T (1 — 2 a) is not balanced anymore (exponential terms have a nonzero 
real part): this rules out the case a = 1/2 and therefore the chance of recovering 
Strichartz estimates with standard arguments from our proof. The analysis of this 
problem will be the object of forthcoming works. 


4. Appendix: special functions, properties and integrals. 

For reader’s convenience we include a short appendix in which we recall defini¬ 
tions and basic properties of the special functions coming into play. We refer to 
[38] and [16] for further details and for a deeper insight of the argument. 

Definition 4.1. Given a±,... ,a p and 61 , 

generalized hypergeometric function as 

p l\{ 1 ■ - - - j Up, b \,..., bq , z ) 


... ,b q complex numbers, we define the 
__ y-"' (u\) n ... (a p ) n z n 

„ (bl)n ■ ■ ■ ( bq) n 77.! 


(4.1) 
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provided the sum of the series is finite, where we are using the Pochhammer symbols 
(a)o = 1 , (a) n = a(a + 1)... (a + n — 1), n > 1 . 


The first thing to investigate is, of course, under which conditions these series 
converge. As a very general framework, we can summarize the situation as follows 

• if p < q the series is convergent for all values of z\ 

• if p = q + 1 the series is convergent for \z\ < 1, and for the special values 
2 = 1 if 


(4.2) 



> 0 


and z — —1 if 



• if p > q + 1 the series never converges except for 2 


0 . 


In this work we are especially interested in the confluent hypergeometric function 
i Fi(a,c,z), which is thus the case p = q = 1. This series is therefore known to 
be convergent for every finite value of z (if c is not zero nor a negative integer). 
Moreover, if a is a negative integer or zero, \F\ reduces to a polynomial of degree 
|a|. Confluent hypergeometric functions can be thought of as a limiting case of the 
celebrated Gauss hypergeometric serie, that namely is case p — 2, q — 1 in (4.1), 
via the relation 


iFi(a,c;z)= lim 2 Fi(a,b]c;z/b). 

6—>-oo 

These functions satisfy the differential equation 


zu" + (c — z)u' — au = 0. 


The asymptotic behaviour of i iq can be derived by writing 

(4.3) i Fi(a,b,z) = ——— - (— z)~ a G(a,a — &+l,z)+ | \ e z z a ~ b G{b — a, 1 — a, z) 


r(b - a ) 

where G has the asymptotic series representation 


T(a) 


G(a, b, z) 


ab 

1+ lb + 


a(a + 1 )b(b + 1) 

21 ? 


+ ... 


Among the properties that confluent hypergeometric functions satisfy we recall 
the following one 


(4.4) 


F(a, c, z) = e z F(c — a, c, —z). 


Among the possible generalizations of hypergeometric function a very useful one 
is given by the so called Appell functions (see [1]). Historically, the motivation 
for introducing series of this type is connected to investigating different kind of 
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products of two different Gauss functions. In particular, we are here interested in 
the following, which is sometimes called Appell function of the second type, 


(4.5) 


oo oo 

F 2 (d; a 1 ,a 2 ;b 1 ]b 2 ]x,y) = ^ ^ 

m= 0n=0 


m ( fl 2)n 

(■ bi) m {b 2 ) n m\n\ 


x m y n . 


It can be seen that convergence condition for F 2 is given by |x| + \y\ < 1. 
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